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Abstract 



We discuss the possibility to extract the neutron-neutron scattering length a nn from experi- 
mental spectra on the reaction 7c? — > iT + nn. The transition operator is calculated to high accuracy 
from chiral perturbation theory. We argue that for properly chosen kinematics, the theoretical 
■ uncertainty of the method can be as low as 0.1 fm. 

o 

^ ■ 1 Introduction 
O 



A precise knowledge of the neutron-neutron scattering length a nn is, e.g., important for an understand- 
ing of the effects of charge symmetry breaking in nucleon-nucleon forces [1]. The scattering length 
a nn characterizes scattering at low energies. It is related to the on-shell 1 Sq scattering amplitude / on 



C$ ' as 

on 1 1 

p r COt5(p r ) - ip r -ann + \r nn pl + 0(pf) - ip r ' 

where p r is the relative momentum between the two neutrons, 5{p r ) the scattering phase shift in the 
1 S'o partial wave and r nn is the effective range. At low energies the terms of order pf. can be neglected 
to very high accuracy. Obviously, a direct determination of a nn in. a scattering experiment is extremely 
difficult due to the absence of a free neutron target. For this reason, the value for &nn is to be obtained 
from analyses of reactions where there are three particles in the final state, e.g. ir~d — > jnn [2, 3, 4] 
or nd — > pnn [5, 6, 7]. There is some spread in the results for a nn obtained by the various groups. In 
particular, two independent analyses of the reaction nd — * pnn give significantly different values for 
a nn , namely a nn = —16.1 ± 0.4 fm [6] and a nn = —18.7 ± 0.6 fm [7], whereas the latest value obtained 
from the reaction ir~d — ► "fnn is a nn = —18.5 db 0.3 fm [4]. At the same time, for the proton-proton 
scattering length, which is directly accessible, a very recent analysis reports a pp = —17.3 ± 0.4 fm 
[8] after correcting for electromagnetic effects. This means that even the sign of A a = a pp — a nn is 
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not fixed. It should be mentioned, however, that state of the art calculations for the binding energy 
difference of tritium and 3 He suggest that A a > [9, 10]. 

In the present work we discuss the possibility to determine a nn from differential cross sections in 
the reaction — > n + nn. Specifically, we show that one can extract the value of a nn reliably by fitting 
the shape of a properly chosen momentum spectrum. In this case the main source of inaccuracies, 
caused by uncertainties in the single-nucleon photoproduction multipole Eq+, is largely suppressed. 
Furthermore there is a suppression of the quasi-free pion production at specific angles. We show that 
at these angular configurations the extraction of a nn can be done with minimal theoretical uncertainty. 

Our investigation is based on the recent work of Ref. [11] in which the transition operator for 
the reaction jd — > ir + nn was calculated up to order x'l 2 in chiral perturbation theory (ChPT) with 
X = m^jM^j ~ 1/7, where m n (Mjy) is the pion (nucleon) mass. Half-integer powers of \ m the 
expansion arise from the unitarity (two- and three-body) cuts (see also [12]). The results of Ref. [11] 
for the total cross section are in very good agreement with the experimental data. The only input 
parameter that entered the calculation was the leading single-nucleon photoproduction multipole Eq + , 
which was fixed from a fourth-order one-loop calculation of Bernard et al. [13]. The uncertainty in 
Eq + is the main theoretical error in the calculation presented in Ref. [11]. Besides this transition 
operator, in the present study we use nucleon-nucleon (NN) wave functions constructed likewise 
in the framework of ChPT, namely those of the NNLO interaction of Ref. [14]. This allows us to 
estimate the theoretical uncertainty which arises from variations in the wave functions. In fact, as 
soon as we include consistently all terms up to order x 5//2 > we expect the ambiguities due to different 
wave functions not to be larger than a % 3 correction, for only at this order the leading counter term 
which absorbs these effects enters. This expectation is indeed quantitatively confirmed in the concrete 
calculations. 

Since we work within chiral perturbation theory we can estimate the effect of higher orders in 
terms of established expansion parameters together with the standard assumption that additional 
short ranged operators, that enter at higher orders, behave in accordance with the power counting 
(the so-called naturalness assumption). This method was also applied in Refs. [15, 16], where the 
reaction 7r~d — > jnn was investigated as a tool to extract a nn . However, to know the effect of higher 
orders for sure, one has to calculate them. Therefore, to derive a reliable uncertainty estimate for the 
extraction of a nn from the jd reaction, we use our leading order calculation as baseline result and 
estimate the theoretical uncertainty from the effects of the higher orders that we calculated completely. 
Based on this, we find a theoretical uncertainty 5a nn < 0.1 fm. We therefore argue that the reaction 
7<i — > ir + nn appears to be a good tool for the extraction of a nn . 

To end this section, we remark that in Ref. [17] a method was proposed to extract scattering 
lengths from jd induced meson production. However, this approach should not be used here, since 
the momentum transfer is not sufficiently large to use this method and with our explicit calculation 
of the transition operator we can reach a significantly higher accuracy. 

2 ChPT calculation for jd — > iz + nn 

The diagrams that contribute to the reaction ^d — ► 7r + nn are shown on Fig. 1. The kinematical 
variables are defined in Fig. 2. 

Before going into the details some comments are necessary regarding the relevant scales of the 
problem. In the near threshold regime of interest here (excess energies of at most 20 MeV above the 

1 Note, that, in contrast to a pp , a nn is not corrected for electromagnetic effects. However, since those are only of the 
order of 0.3 fm [1] they are not relevant for the sign of A n . But they ought to be taken into account for determining 
charge symmetry breaking effects quantitatively. 
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Figure 1: Diagrams for 7c? — * 7r + nn. Shown are one-body terms (diagram (a) and (b) ), as well as the 
corresponding rescattering contribution (c) — all without and with final state interaction. Diagrams 
(d) shows the class of diagrams with intermediate NN interaction. Solid, wavy, and dashed lines 
denote nucleons, photons and pions, in order. Filled squares and ellipses stand for the various vertices 
(see Ref. [11] for the details), the hatched area shows the deuteron wave function and the filled circle 
denotes the nn scattering amplitude. Crossed terms (where the external lines are interchanged) are 
not shown explicitly. 



pion production threshold) the outgoing pion momenta are small compared even to the pion mass. 
Thus, in addition to the conventional expansion parameters of ChPT m n /A x and q^/A x , where A x 
denotes the chiral symmetry breaking scale of order of (and often identified with) the nucleon mass, 
and q 7 denotes the photon momentum in the center-of-mass system which is of order of the pion 
mass, we can also regard k^/m^ as small, where k n denotes the momentum of the outgoing pion. In 
what follows we will perform an expansion in two parameters, namely 

Xm = m w /M N and xq = K/m n . 

Obviously, the value of the second parameter depends on the excess energy Q. The energy regime 
of interest to us corresponds to excess energies up to 20 MeV. The maximum value of %Qi XQ ax = 
^2Q/m- K , at the highest energy considered is thus about 1/2. Since this is numerically close to ^JXm 
we use the following assignment for the expansion parameter: 

X ~ Xm ~ Xq • (2) 

The tree level 7p — * Tr + n vertex, as it appears in diagrams (al) and (a2) in Fig 1 (the vertex is labeled 
as filled square), contributes at leading order (order x )) and orders x 1 and X 2 -, depending on the 
one-body operator used. Note that the loop diagrams with irN rescattering (see diagrams (b), (c) 
and (d) in Fig. 1) contribute at order x m as well as at x m XQi Xm, 2 and at xII^Xq- The origin of the 
non-integer power of x are the two-body (7riV) and three-body (nNN) singularities. Thus, all terms 
up to x 5 ^ 2 are explicitly taken into account in our calculation of the transition operator. 

As already emphasized, we employ wave functions evaluated in the same framework in order to 
have a fully consistent calculation. In our work, we use the N 2 LO wave functions corresponding to 
the chiral NN forces introduced in Ref. [18] and based on the spectral function regularization (SFR) 
scheme [19]. At this order, the NN force receives contributions from one-pion exchange, two-pion 



3 



d ^ P2 n 

Figure 2: Kinematical variables for 7c? — > ir + nn. The relative neutron-neutron momentum is denned 
as p r = \(pi -p 2 ). 



exchange at the subleading order as well as from all possible short-range contact interactions with up 
to two derivatives. In addition, the dominant isospin-breaking correction due to the charged-to-neutral 
pion mass difference in the one-pion exchange potential together with the two leading isospin-breaking 
S-wave contact interactions were taken into account [18]. The two corresponding low-energy constants 
were adjusted to reproduce the scattering lengths a nn and a pp . The SFR cutoff A is varied in the range 
500 . . . 700 MeV. It was argued in Ref. [19] that such a choice for A provides a natural separation of 
the long- and short-range parts of the nuclear force and allows to improve the convergence of the chiral 
expansion [19]. The cutoff A in the Lippmann-Schwinger equation is varied in the range 450. . . 600 
MeV. For an extensive discussion on the choice of A and A the reader is referred to [14, 18]. 



3 Differential cross sections: relevant features 

In this section we outline the features of the differential cross section for unpolarized particles that are 
important for our considerations. For later convenience let us consider the function F proportional to 
the square of the matrix element as well as the five-fold differential cross section 

d 5 a(p r ,9 r ,(j) r ,9 n ,(p n ) 



F(p r ,9 r ,(j) r ,6 n ,(l> n ) = Cp r fc,r(Pr)|A1(p r ,0 r ,0 r ,0,r,0,r)| 2 K jj=fZ dQ^ dpi ' ^ 

where p r (k w ) stands for the relative momentum of the two final neutrons (momentum of the final 
pion) in the center-of-mass frame, r , cp r (9 n , 4> n ) for the corresponding polar and azimuthal angles, 
respectively, and \M\ 2 for the squared and averaged amplitude. In Eq. (3) C is an irrelevant dimen- 
sioned constant. In what follows we will consider only shapes of cross sections and therefore the value 
of C is not important for our considerations. The value of k w at given p r and excess energy Q is fixed 
by energy conservation: 

„2 L.2 u2 

O = — + n + — — (4) 

v M N AM N 2m w ' V ' 

hence we write ^(pr) in Eq. (3). 

In the following we choose the momentum <f 7 of the initial photon to be along the z-axis. Then 
the cross sections at a certain excess energy Q depend on four variables, namely the magnitude of 
the relative momentum of the two final neutrons p r , the polar angles of the vectors p r and k n , and 
the difference between the azimuthal angles of those two momenta. Unpolarized cross sections are 
invariant under rotations around the beam axis, which makes the dependence on the missing angle 
trivial. 
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Typical differential cross sections F are shown in Fig. 3 as a function of p r at some fixed set of 
angles {4> r , 9 n , 4> w } and Q = 5 MeV for two different values of 9 r . One can see from this figure that 
for the differential cross section F of Eq. (3) there are two characteristic regions: 

1. The region of quasi- free production (QF) at large p r , which corresponds to the dominance of those 
diagrams of Fig. 1 that do not contain the NN interaction in the final or intermediate states. 
In the Appendix we give explicit expressions for the diagram al - the most significant diagram 
of this type. At large p r the pion momentum k n is small (see Eq. (4)) and the arguments of the 
deuteron wave function in Eqs. (A.l) and (A. 2) may become small for particular combinations 
of ±p r and q-,/2. This feature gives rise to a peak in the differential cross section at large p r . 

2. The region with prominence of the strong nn final-state interaction (FSI) at small p r (in fact, 
we would have the strongest final state interaction at zero relative momentum, however the cross 
section goes to zero at p r = due to the phase space, therefore we see a peak shape). 

One can see from Fig. 3 that the FSI peak depends on the value of 9 r only marginally, whereas the 
quasi-free peak shows significant dependence on this angle. In particular, the quasi-free production 
is largely suppressed at 9 r = 90° — at this angle the arguments of the wave functions in both terms 
in the r.h.s. of Eqs. (A.l) and (A. 2) are large. It can also be seen from Fig. 3 (right panel) that 
the effect of higher orders is more important for the quasi-free production amplitude — the influence 
of higher-order effects on the FSI production is quite small. Another interesting observation is that 
the contributions of higher orders change the relative height of the two peaks - the FSI peak goes up 
whereas the QF peak goes down when we proceed from the LO calculation to the order X > / 2 . In order 
to suppress the distortions of the spectrum due to higher orders in the chiral expansion, which is the 
condition for an extraction of a nn with small theoretical uncertainty, configurations should be chosen 
where 9 r = 90°. 

We now briefly discuss the dependence of the cross section on the remaining angles 9 W , <j) w (we always 
may choose <j> r to be zero). The dependence on 9 n is illustrated in Fig. 4. One can see from this figure 
that the dependence on 9 n is significant for both the quasi-free as well as the FSI peak. This can be 
easily understood from the explicit expressions for the matrix elements given in the Appendix keeping 
in mind that already at Q = 5 MeV the maximal value of k n is about m n /3 while q 1 « m n . Thus, 
the momentum transfer to the nucleon pair, \q 7 — k n \, varies in the range 2m 7r /3 to 4m 7r /3 depending 
on 9 n . Since the <S-wave deuteron wave function is large only for very small arguments, the influence 
of the direction of k n is significant. In addition, from Fig. 4 it follows that a variation of 9 n not only 
changes the magnitude but also the shape of the cross section, even in the FSI region. This has to be 
taken into account in the analysis of any experiment. 

In contrast to the polar angles, the dependence of F on 4> n is negligible for all configurations (there 
is no dependence at all for 9 r = 0° and at 9 r = 90°, only the anyway small QF contribution changes 
by just 5 %). 

4 Extraction of a nn and estimate of the theoretical uncertainty 

In this section we discuss how to extract the scattering length from future data on 7c? — > ir + nn as well 
as the resulting theoretical uncertainty. Our focus is especially the latter point. As in the previous 
section we will only discuss results at excess energy Q = 5 MeV. However, the analysis can be repeated 
analogously at any excess energy within the range of applicability of the formalism, i.e. Q < 20 MeV. 

We are interested in extracting the value of a nn , which, in turn, is a low-energy characteristic of 
neutron-neutron scattering and manifests itself in the momentum dependence of the cross section at 
small values of the momentum p r . The influence of the value of a nn on the cross section is illustrated in 
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Figure 3: Left panel: Differential cross section. The solid line corresponds to the configuration 
when the quasi-free peak is suppressed (0 r = 90°), whereas the dashed line corresponds to one of 
the configurations when the quasi-free production amplitude is maximal (6 r = 0°). The values of 
the remaining angles are 6^ = 135°, <p r = 4> n = 0°; they are the same for both curves. Right panel: 
Differential cross section — relative strength of QF and FSI peaks. Here the dashed curves correspond 
to the calculation at LO, the solid ones to the calculation at x 5/2 - Curves denoted by "FSI" ("QF") 
are obtained by retaining only those diagrams of Fig. 1 that contain (do not contain) the final or 
the intermediate nucleon-nucleon interaction. The labels "0 degrees" and "90 degrees" denote the 
corresponding values of 6 r for the "QF" curves whereas the "FSI" curves are almost insensitive to this 
angle. The values or the remaining angles are as on the left panel of this figure. The overall scale is 
arbitrary in both panels but the relative normalization is the same for all curves. 

Fig. 5, where the cross sections are shown for three different values of a nn , namely —18, —19, — 20 fm. 
For each value there are two curves, the dashed one corresponds to 9 r = 0°, and the solid one to 
6 r = 90°. One can see from Fig. 5 that the influence of different values of a nn is significant in the FSI 
peak and marginal in the quasi-free peak, as one would have expected. 

In the previous section we have shown (see right panel in Fig. 3) that the relative height of the quasi- 
free and the FSI peak changes if the effects of higher orders are included in the cross sections. Therefore 
those angular configurations are to be preferred, where the quasi-free production is suppressed. 

The central point of this study is to demonstrate that there is a large sensitivity of the momentum 
spectra to the scattering length and that this scattering length can be extracted with a small and 
controlled theoretical uncertainty. As outlined in the Introduction, we can estimate this uncertainty 
reliably, because the effect of the higher orders up to x 5 ^ 2 are calculated completely. In order to 
demonstrate the effect of those higher orders on the shape of the momentum distribution, in Fig. 6 
we show as the light band the spread in the results for the calculation from LO to x 5 ^ 2 - The results 
also include higher partial waves for the pion as well as the final nn system. 

There is some sensitivity to the behavior of the deuteron wave function at short distances. For the 
reaction ir~d — * jnn this sensitivity was identified as the largest effect at N 3 LO in Ref. [15] 2 . Guided 
by that observation we include in the uncertainty estimate also the spread in the results due to the 
use of different wave functions. In order to remove the effect of the change in normalization when, 
e.g., changing the chiral order, all curves are normalized at p n orm = 30 MeV in Fig. 6. In the same 

2 Within the framework of ChPT with a consistent power counting scheme, the quantitative impact of the wave- 
function dependence is governed by the order at which a counter term appears that can absorb this model dependence. 
The corresponding counter term for the *yd as well as the ivd reaction arises at N 3 LO. 
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Figure 4: Dependence of the differential cross section on 8 n . The left panel corresponds to the 
suppressed quasi-free amplitude (6 r = 90°), the right panel to the maximal quasi-free amplitude 
(9 r = 0°). Solid, dashed, dotted, and dash-dotted lines correspond to 6 n = 0°, 45°, 90°, 135° 
respectively. The values of the remaining angles are <p r = 4> n = 0°. The overall scale is arbitrary in 
both panels but the relative normalization is the same for all curves. 

Figure (with the same normalization) we also show the change in the shape that comes from different 
values of the scattering length: the dark band is generated by a variation of a nn by ±1 fm around 
the central value of — 18.9fm. Clearly, the theoretical uncertainty is negligibly small compared to the 
signal of interest. 

One way to quantify the theoretical uncertainty is through the use of the function S, defined as 

Pmax 

S(a nn ,$) = J dp r (V(p r |a^,$ (0) ) - N(a nn ,<S>) F(p r \a nn ,$)) 2 w(p r ) , (5) 
o 

where p max = \/MjyQ is the maximum value of p r , F(p r \a nn ,<&) is proportional to the five-fold 
differential cross section as defined in Eq. (3). In the latter we refrained from showing the angular 
dependence in favor of the parametric dependence of the cross section on the titi scattering length d nn 
as well as the multi-index <£, which symbolizes the dependence of the cross section on the chosen chiral 
order and the wave functions used, as outlined above. The weight function w{p r ) was introduced to 
allow us to suppress particular regions of momenta in the analysis — the role of w(p r ) will be discussed 
in detail below. For simplicity we may assume that S is dimensionless; all dimensions can be absorbed 
into the constant C defined in Eq. (3). 

The value ann denotes the central value of the scattering length (—18.9 fm) for which we perform 
the estimate of the theoretical uncertainty 3 whereas <3?(°) corresponds to the baseline type of calcu- 
lation, namely leading order with chiral wave functions as specified in the Appendix. The relative 
normalization N(a nn , $) is fixed by demanding that S gets minimized for any given pair of parameters 
a nn ,<& (dS/dN = 0). This gives 

Pmax , * 

/ dp r F(p r \a n n,$ i0) )F(p r \a nn , $)w(p r ) 
N ( a nn, *) = -2 ^ x • (6) 

/ dp r F 2 (p r \a nn ,$)w(p r ) 




^Note that the theoretical uncertainty practically does not change when the central value of the scattering length 
varies in the relevant interval ±1 fm. 
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Figure 5: The effect of varying the value of a nn on the differential cross section. The solid and dashed 
lines correspond to the same angular configurations as in Fig. 3, left panel. The different values of 
a nn are shown on the figure. The overall scale is arbitrary but the relative normalization is the same 
for all curves. 



Obviously S is the continuum version of the standard x 2 sum, i.e. it characterizes the mean-square 
deviation from the baseline cross section F(p r \a^,^). In this way we determine the theoretical 
uncertainty in full analogy to the standard method of data analysis. 

In order to quantify the theoretical uncertainty we may define $ max as that chiral order and choice 
of wave function, where S(a,nn, ^max) gets maximal: 

<D max ) = max {s(o£L *)} • (7) 

Therefore 5(o™,$ max ) provides an integral measure of the theoretical uncertainty of the differential 
cross section. Demanding that the effect of a change in the scattering length by the amount Aa nn 
matches that by the inclusion of higher orders etc., we can identify Aa nn as an uncertainty in the 
scattering length. Expressed in terms of S, we may define Aa nn via 

S(a<S> + Aa nn , $(°>) = S(a&), <& max ) . (8) 

This relation is illustrated in Fig. 7. The dashed horizontal line corresponds to 5(c4m, $ mm ), where 
we use w{p r ) = 1. The dashed parabolic line shows the corresponding S{anl + Aa nn , <I>( )) as a 
function of Aa nn . The calculation is performed for 8 r = 90°, and 6^ = 0°. The crossing point of the 
curves corresponds to Aa nn = 0.16 fm, which can be identified as the theoretical uncertainty for the 
extraction of the scattering length. 

In the previous section we showed that the signal region is located at momenta lower than 30 MeV. 
On the other hand, the theoretical uncertainty of the differential cross section is largest for large values 
of p r due to the onset of the quasi-free contribution. In view of these two facts it seems reasonable 
to use such weight functions w(p r ) that suppress the contribution of large momenta. For instance, we 
may use w(p r ) = @(p cut — p r ) for the weight function. If we choose, e.g., p cut = 30 MeV the theoretical 
uncertainty of the extraction of the scattering length reduces to 0.07 fm, as is demonstrated by the 
solid lines in Fig. 7. This figure nicely illustrates that the parabolic curve that represents the signal 
changes only very little when a restriction to small values of p r is applied. At the same time this 
procedure significantly reduces the value of the uncertainty S(c4°2, •I'max)- 
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Figure 6: The light (white) band is the error band, and dark (blue) band correspond to ±1 fm shift 
in the scattering length from the central value —18.9 fm. 

The observation that the dependence of 5(ai°2 + Aa nn , <1>( )) on Aa nn is very well approximated by 
a parabola allows for a more systematic study of the p cut dependence of the theoretical uncertainty. 
We therefore define 

_ S(a^ + Aa nn ,^)\ P ^) 
[P j " (Aa™) 2 ' W 

where the explicit p cut dependence is introduced into the function S through the weight function w as 
explained above. The dashed and the solid parabola in Fig. 7 can then be written as a(p cnt ) (Aa nn ) 2 , 
with a(p max ) = 41 fm -2 and a(30 MeV) = 33 fm -2 . In the left panel of Fig. 8 we show a(p cnt ) as 
the solid line. In the same panel the dashed line represents the measure of the theoretical uncertainty 
given by 5(al°2, <3? max |p cut ), multiplied by a factor of 40. This figure makes more quantitative the 
statement made above: for very small values of p cut we cut into the signal region and therefore a 
shows a very rapid variation. However, as soon as p cut is larger than 30 MeV it goes to a plateau (in 
the figure indicated by the arrow). On the other hand, the theoretical uncertainty is monotonously 
growing once p cut is larger than 30 MeV. From this figure we deduce that the ideal value for p cut 
is between 25 and 40 MeV. This translates into a theoretical uncertainty between 0.05 and 0.1 fm, 
as illustrated in the right panel of the same figure. The value of 9^ also has some impact on the 
theoretical uncertainty, however, in its whole parameter range the estimated uncertainty stays below 
0.1 fm for p cut = 30 MeV. 

Clearly, also the experimental data, once they exist, should be analyzed using a procedure analogous 
to the one given above. This means that the scattering length is to be extracted from a x 2 fit of the 
theoretical curves to the data. In this work we used the calculation at LO as baseline result and the 
results at higher orders to estimate the theoretical uncertainty. Consequently, we propose to use the 
momentum spectrum calculated at LO in the fitting procedure of the experiment. The corresponding 
analytical expressions are given in the Appendix. The only parameter to be adjusted besides the 
scattering length is the overall normalization. In this fitting procedure only those data points should 
be included that are below a given p cut , in order to keep the theoretical uncertainty small. 
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Figure 7: The functions 5(a„,i, <5 max ) and S(dnn + Aa nn ,$^) are shown by the horizontal and 
parabolic curves, respectively. The solid curves are obtained by adding the weight factor in Eq. (5) 
that cuts all momenta above 30 MeV in distinction from the dashed ones. The calculation is performed 
for the scattering length a,nn = —18.9 fm, 9 r = 90°, and 9^ = 0°. The value of Aa nn corresponding 
to the crossing point of the horizontal and parabolic curves determines the theoretical uncertainty of 
the calculation. 



5 Discussion and conclusions 

In the previous section it was shown that for the angular configurations that suppress the quasi-free 
production the inclusion of higher order effects (NLO, N 2 LO, and x 5 ^ 2 ) as wen as the use of different 
wave functions leads only to a minor change in the momentum dependence of the five-fold differential 
cross sections. 

Based on this observation we propose to use the momentum spectrum calculated at LO for the 
extraction of the neutron-neutron scattering length from the data. This procedure has the advantage 
that the corresponding matrix elements can be given in an analytic form (see Appendix) that could 
be used directly in the Monte Carlo codes for the experiment analysis. In this way the non-trivial 
dependence of the spectra on 6^, discussed above, can be easily controlled. The scattering length can 
then be extracted by a two parameter fit to the data where, simultaneously to a variation in a nn , the 
normalization constant needs to be adjusted. 

Note that the leading order calculation basically agrees to the expression given in Ref. [20] long 
ago. However, a systematic and reliable study of the theoretical uncertainties of the extraction was 
possible only within our full calculation up to order x 5 / 2 in ChPT. In this way we could show that 
the reaction jd — ► Tr + nn is very well suited for a determination of the nn scattering length. The 
theoretical uncertainty of order 0.1 fm for the extracted scattering length, estimated in this paper, is 
of the same order as that claimed for ir~d — > r ynn [16] and nd — > pnn [4, 7]. 

We discussed in detail the theoretical uncertainty for a fixed excess energy of Q = 5 MeV only, 
however, it should be clear that the procedure can be easily repeated for any energy within the range of 
applicability of our approach [Q < 20 MeV). For example, we checked that the theoretical uncertainty 
stays below 0.1 fm also at Q = 10 MeV. Note that the number of events in the signal region scales 
roughly with y/Q, the phase space available for the pion. It remains to be seen which energy is the 
best for the corresponding experiment. 

We showed that for a proper choice of both kinematics and weight function uu, the theoretical 
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Fi gure 8: Left panel: Comparison of the p cn ^ dependence of functions S(cinn, ^max 

\p cut ) (dashed curve) 

and a(p cnt ) (solid curve). The calculation is performed for the scattering length cinn — -18.9 fm, 
9 r = 90°, and 9 n = 0°. Right panel: The corresponding theoretical uncertainty Aa nn as a function of 

pCut 



uncertainty for the extraction of the neutron-neutron scattering length from jd — ► TT + nn can be as 
low as 0.1 fm. It should be stressed, however, that this error was evaluated most conservatively - we 
use our LO calculation as baseline result and estimate the theoretical uncertainty from the effects of 
the higher orders that we calculated completely. This error can be significantly reduced by further 
studies. For example, if we include in the uncertainty estimate only the spread in the results due to 
the use of different wave functions, which is identified as the largest effect at N 3 LO for the reaction 
ir~d — > jnn [15], the theoretical uncertainty of the extracted scattering length reduces by one order 
of magnitude. This indicates that the theoretical uncertainty is indeed under control. However to 
put this N 3 LO estimation on more solid ground a complete calculation should be performed to this 
order. Most of the operators that are relevant at this order are the same as those of ir~d — ► jnn, given 
explicitly in Ref. [21]. One counter term enters, which can be fixed from other processes [16], e.g., 
from nd scattering [22], the reaction NN — * NNir [23], or from weak decays [16]. Once this is done 
we may use our calculation to order x 5 ^ 2 as baseline result and estimate the theoretical uncertainty 
from the then available N 3 LO calculation. 

Although we have identified the angles 9 r = 90° as the preferred kinematics, also other configura- 
tions could be studied in order to control the systematics. However, then the spectra calculated at 
X 5//2 should be used in the analysis. 
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A Leading amplitudes 



In this appendix we give explicit expressions for the amplitudes that appear at leading order in the 
calculation for — > i: + nn. As outlined in the main text these expressions can be used directly in 
the analysis of the data, once available. In addition, they should also proof useful for the design of 
the experiment. Note, as outlined in the text, only near 9 r = 90° the leading order calculation gives 
a sufficiently accurate representation of the spectra. At all other angles one should use the complete 
calculation. 

At leading order only diagrams al and a2 of Fig. 1 contribute. Since only the momentum dependence 
of the amplitudes is relevant for the experimental analysis we drop an overall factor compared to 
Ref. [11]. The corresponding amplitudes read: 

Ki = - 4/2 + <f 7 /2) + u{-p r - 4/2 + g 7 /2)) (A.l) 

Ki = ( u (Pr ~ 4/2 + ? 7 /2) - u(-p r - K/2 + <f 7 /2)) (A.2) 

w a f° n (Pr) f d 3 p u(p-k w /2 + <f 7 /2) g(p) 

M a2 = 87T 



g{Pr) J (2vr) 3 p 2 -p 2 .- iQ 

f ml {Pr) C i D j ln ( a i ~ { Pr + i%n _ + Pj ^jgvA ^ ^ 

where u(p) denotes the 5-wave part of the deuteron wave function in momentum space. We checked 
by explicit calculations that the inclusion of the deuteron D-w&ve changes only the absolute scale of 
the differential cross sections but not its momentum dependence. Thus, the D-wave contribution is 
not taken into account in the parameterization. The quantity is defined as q nl = \k n — <f 7 |/2. The 
labels s and t stand for spin singlet and triplet final two-nucleon states, respectively — we do not 
write out the corresponding spin structures. We take into account only the 1 So partial wave in the 
final state interaction. For a discussion of the effect of nn P-waves see Ref. [11]. 

To derive the expression for M a2 we used the fact that the neutron-neutron scattering amplitude can 
be represented to high accuracy in separable form [11, 24]. The neutron-neutron scattering amplitude, 
f(p, k; E), can be written in half off-shell kinematics as 

f(p,k-y,M N)= ^^f'f', =r(kM, (a.4) 

where the corresponding on-shell amplitude f on (k) can then be expressed in terms of the scattering 
phase-shifts through 

For small momenta one can use the effective range expansion for k cot 5 = —l/a nn + r nn k 2 /2 + 0(k 4 ), 
in agreement with Eq. (1). Here a nn is the parameter to be fitted to the data and r nn = 2.76 fm. We 
checked that changing the value of r nn within the bounds allowed (±0.1 fm [1]) leads to negligible 
effects on the extraction of the scattering length. In this way we expressed the matrix element explicitly 
in terms of the scattering length. We checked that the ratio g{p)/g{k) in Eq. (A.4) does not change 
when we vary the scattering length within acceptable range bounds. 

In order to evaluate the convolution of the deuteron wave function with the nn final state interaction 
analytically, we needed to employ the following parameterizations for the 1 So nn form factor g(p) (see 
Eq. (A.4)) and the S-wave deuteron wave function 



i 
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l Sq form factor 


S-wave deuteron w.f. 




Pi [MeV] 


A [MeV] 


oii [MeV] 


Ct [MeV 1 / 2 ] 


1 


164.53278 


31.101228 


45.334919 


43.543212 


2 


246.85751 


-1310.3056 


242.66091 


-35.643003 


3 


329.18224 


9455.9603 


439.98691 


419.25214 


4 


411.50697 


-9666.0268 


637.31291 


-1833.4708 


5 


493.83170 


-55571.615 


834.63891 


-3710.8173 


6 


576.15643 


64600.071 


1031.9649 


24903.150 


7 


658.48116 


149128.85 


1229.2909 


-31673.576 


8 


740.80589 


-84844.967 


1426.6169 


26476.636 


9 


823.13062 


-295594.17 


1623.9429 


-118733.48 


10 


905.45536 


-30332.710 


1821.2689 


259759.15 


11 


987.78009 


560829.89 


2018.5949 


-223816.07 


12 


1070.1048 


-307006.25 


2215.9209 





Table 1: Parameters of the So form factor and the S-wave deuteron wave function for the separable 
representation of the N 2 LO chiral NN potential. 



where the parameters corresponding to the ChPT calculation at N 2 LO with cut offs {A, A} = {550 MeV, 
600 MeV} (see Ref [18] for details) are listed in Table 1. Note that the coefficients in the parameter- 
ization of the wave function have to fulfill the relation ^ Cj = in order to ensure the regularity of 
the deuteron wave function at the origin in coordinate space [25]. 

The squared and averaged amplitude to be used in the expression for the differential cross section, 
defined in Eq. (3) is 

\M( Pr ,e r , </> r , e n , (Mfitl 2 = \Ki + M a2 \ 2 + 2 \M l al | 2 . (A.e) 

In a fit to data two parameters are to be adjusted, namely the overall normalization C of Eq. (3) 
and the object of desire, a nn . 
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